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WEIGHTED POLYNOMIAL APPROXIMATION 
OF RATIONAL BEZIER CURVES* 


STANISEAW LEWANOWICZt, PAWEL WOZNY, PAWEL KELLER 


Abstract. We present an efficient method to solve the problem of the constrained least 
squares approximation of the rational Bezier curve by the Bezier curve. The presented al¬ 
gorithm uses the dual constrained Bernstein basis polynomials, associated with the Jacobi 
scalar product, and exploits their recursive properties. Examples are given, showing the 
effectiveness of the algorithm. 

AMS classification-. Primary 41A10. Secondary 65D17, 33D45 

Keywords: Rational Bezier curve; Polynomial approximation; Constrained dual Bernstein 
basis. 


1. Introduction 


In CAGD, it is frequently important to approximate a rational Bezier curve by a poly¬ 
nomial one. In the last two decades, many approaches to this problem have been proposed 
m na m 0 m im n\. The large spectrum of methods contains approximation by hy¬ 
brid curves [T8, IS [15], Hermite interpolation [5], [15], progressive iteration approximation 
HU, least squares approximation [10] and approximation by Bezier curves with the control 
points obtained by successive degree elevation of the rational Bezier curve [7]. 

In this paper, we consider the following approximation problem. 

Problem 1.1. Let be given a rational Bezier curve R n of degree n, with control points 
ri G and positive weights u>i G M, 




where 


B?(t) = f( 1 - (0 < i < n) 

are Bernstein basis polynomials. Find a degree m Bezier curve 



m 


( 1 . 2 ) 



(0 < t < 1) 


i =0 
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such that the value of the error 

(1.3) f\l-t) a t^\\R n (t)-P m (t)\\ 2 dt (a, (3 > —1) 

Jo 


is minimized in the space TYf of parametric polynomials in W l of degree at most m (for 
simplicity, we write H m := U l m ) under the additional conditions that 


(1.4) 


Pi i} (0) = P®(0) (i = 0,l,...,fc-l), 
pW(l) = F^(l) (j = 0,1,..., Z — 1), 


where k + / < m. Here 


is the Euclidean vector norm. 


Note that in the case a = (3 = 0, the above problem as well the method proposed in this 
paper reduce to the form given in jl4| . 

The values of n and m are not related. However, if (m + 1 )d < {n + l)(d + 1) then 
the number of parameters of the approximating polynomial curve is smaller than the total 
number of parameters of the rational curve. 

In |TU] . an approximate solution to the above problem with m > n is obtained by solving 
a linear least squares problem. Paper [2] deals with more general problem of the constrained 
degree reduction of rational Bezier curves; one of the auxiliary problems discussed there 
contains the above problem as a particular case. We present a method which is based on 
the idea of using constrained dual Bernstein polynomial basis [20] to compute the control 
points pi . Our algorithm is efficient thanks to using fast schemes of 1° evaluation the 
Bezier form coefficients of the dual polynomials [13] and 2° numerical computation of the 
collection of integrals 


Jo UJ V') 

The cost is significantly lower than in the case of the special variant of the method of [2], 
which needs inverting a (n + m — k — Z + l) x (n + m — k — Z + 1) matrix. 

Let us mention that the problem stated above can be also considered for other norms. 
However, even the simpler problem of constrained degree reduction of Bezier curves in 
Loo-norm requires higher computational complexity (see, e.g., &)■ The most appropriate 
metric for curves in geometric terms would be the Hausdorff distance, but the computation 
of such distance of the nonlinear curves is not so easy. Hence, I/ 2 -norm seems to be a good 
choice as we can construct solution in explicit form using the Bernstein and dual Bernstein 
bases in a natural and convenient way. 

The outline of this paper is as follows. Section [2] contains basic facts on the constrained 
dual Bernstein polynomials. Section [3] brings a complete solution to Problem 11.11 for 
implementation details, see Section [4] In Section [5], the proposed method is applied to 
some examples and compared with two other algorithms. 
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2. Constrained dual Bernstein polynomials 

Let where k + l < m, be the space of all polynomials of degree < m, whose 

derivatives of order < fc — 1 at i = 0, as well as derivatives of order < l — 1 at t — 1, vanish: 

nW) := {P e n m : P«(0) = 0 (0 < i < k~ 1) andP (j) (l) - 0 (0 < j < l - 1)} . 

Obviously, dim II,= m — k — l + 1, and the Bernstein polynomials {B™, B ™ +1 ,..., 
form a basis of this space. There is a unique dual constrained Bernstein basis of degree m, 

D { ™’ k ' l) (x- a, 0), D%$' l) (x; a, /3),..., D ( ^ k l ,l \x- a, 0), 

satisfying 

B T) = S ‘i = 

where 8^ is 1 if i — j and 0 otherwise, and the inner product (•, •) is given by 

(/, g) ■= [ (i - dt {a, (3 > -1). 

Jo 

Properties of the polynomials are studied in [T2J [T3j and [20]. We need the 

following result. 

Theorem 2.1 ([T3]). The constrained dual basis polynomials have the Bezier-Bernstein 
representation 

m—l 

(2.1) D ( f n ’ k ’ l \x] a, (3) = Cij(m, k, l, a, (3) B™(x), 

j=k 

where the coefficients Cij = Cij(m,k,l,a, j3) satisfy the recurrence relation 

Ci+ij =yr 7 W {(i - j)(2i + 2j - 2m - a + (3) a, 

A \ l ) 

+B(j) c i} j -1 + A(j) Qj+i B{i) Ci-ij} 

(k < i < m — l — 1, k < j < m — l) 

A{u) := (u — m)(u — k + l)(n + k + f3 + l)/(u + 1), 


( 2 . 2 ) 


with 


B{u ) := u{u — m — l — a — l)(w — m + l — 1 )/(u — m — l). 

We adopt the convention that c^ := 0 if i < k, or i > m — l, or j < k, or j > m — l. The 
starting values are 

'rn\ (—l) fc (<r + 2k + 21 + 1 ) m -k~i(k + (3 + 2 ) m -i 
k ) (m — k — l)\B{a + 21 + 1, (3 + 2k + 1) 


(2.3) 


c kj 


fm — k — l\fm\ 1 (“l )' 7 

\ j ~ k ) \ j J {a + 21 + l) m -i -j(k + (3 + 2)j 





4 


STANISLAW LEWANOWICZ f , PAWEL WOZNY, PAWEL KELLER 


where j = k,k + 1,... ,m — l, and we use the notation 
k -1 

{a)k := JJ(a + j) (k > 0), cr := a +/3 + 1, 
j=o 


B( A, //) 


r(A)rfc) 
r(A + /i) 


Observe that the quantities c^- can be put in a square table (see Table |T|). Now, the 


Table 1. The c-table 


0 0 

0 Cfcfc Cfc,fc-1-1 

0 Cfc+i,fc ,fc+i 


0 

Ck,m—l 0 
Qc4T,m —l 0 


0 Cm—Z,fc 


0 


0 


Cm— l.rri —/ 0 


0 


c-table can be completed very easily in the following way. 

Algorithm 2.2 (Computing the coefficients Cij(m,k,l,a, /3) [T3]). 

I: Compute recursively quantities Ckk, Cfc,fc+i, ■ ■ ■, Ck :m -i, filling the first row of the 
c-table, by the formulas 

+ 2k + 2l + !)„*, 


(2.4) 

(2.5) 


Cfc m—Z • 


1 / \-l 

m\ m 


ij B(a T 21 T 1, /3 T 2k + 1) (m — k — l )\’ 


c kj ■ 


(.3 ~ m){j - k + l){j + fi + k + 2) 




(j + l)(j - "i + l)(j - a - l - m) 

(j — m — l — 1 , m — l — 2,.. ., k). 

II: For i = k, k + 1,... ,m — l — 1 and j = k,k + 1,... ,m — l, compute c i+1 j, using 
the recurrence (12.21) . 

3. Constrained polynomial approximation of rational Bezier curve 


Clearly, the Bezier curve being the solution of Problem 11.11 can be obtained in a com¬ 
ponentwise way. Hence, it is sufficient to give the details of our method of solving this 
problem in case where R n , P m G R 1 . 

Given the rational function 

Qn(t) _ Er=o 


Rn(t) = 


u(t) E"=0 Vi B i{t) 


(3.1) 


(0 < t < 1), 
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with fj 6 K and E R+, we look for a degree m polynomial 

m 

(3.2) P m (t) = ^r,piB™(t) (0 < t < 1) 

i=0 


which gives the minimum value of 

(3.3) ll-Rn(^) - Pm(t)\\ 2 L 2 = (Rn ~ Pm, Rn ~ Pm) 
with the constraints 

id°(0) = Pii ) (0) (i = 0,1,..A: — 1), 

(3.4) 

R^\l) = P^\l) (j = 0,1,..., / — 1), 

where k + l < m. 


Theorem 3.1. Given the coefficients r 0 ,r 1 ,... ,r n and weights co 0 , or,..., u) n of the rational 
function (13.11) . the coefficients Po?Pi? • • • ,p m of the polynomial (13.2)1 minimising the error 
with constraints (13.41) are given by 


(3.5) 


Pi 


(m — i)\ 


i—l 


m\ 


3=0 


Sm ~ Afi (s = 0,1,..., 1); 


(3.6) Pm _,=(-i (i = 0,l,...,(-l); 


mi 


3 = 1 


(3.7) 


Pi = 


EdK- E 


m—l / \ / \—1 

m \ ( n + m 


j J \ J + h 


Cij(m,k,l,a,/3)I j+h 


( k —1 m \ 

S ) P.i K ii (i = k,k + l,...,m-l), 

j =0 j=m—l+LJ 

where Cij(m,k,l,at, j3) are introduced in (12.11) . while and o t \ are defined recursively for 
i = 0,1,... by 


(3.8) 


n\ 

Qio '■ = — 
u 0 


i —1 


1 Ai/ A / M ^ 3UJ ° 

7-ttA (cu 0 r 0 ) - > . 7 -r—rr fto 

(«-*)! pj W (n-i+j)\ 


Qil • 


n! 


CJrj 


i —1 


(n — *)! 






pj W (^-* + j) ! 


I ^J ' 1 


(3.9) 
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(3.10) 

(3.11) 


h , = f 1 (i - dt 

h ' Jo w(t) 

K =( m \( m V 1 ~ j)m-k-l +1 

lJ V?/ W {i - ]){i - k)\{m - l - i)\ 


(a + l + l) m _j(/3 + k + l)j 
(a + / + l) m _i(/3 + k + l)j 


//ere we use the standard notation A°c/ l = q>. A j c/j = A J 1 c/ l+ i — A J 1 c/ l (j = 1,2,.. .j. 


Proof. Recall that for arbitrary polynomial of degree N , 

N 

u N (t) = j2 UiB ?( *)> 

i=0 

the well-known formulas hold (see, e.g., d p. 49]) 


,.v AT! ATI 

£$ } (0) = •„. A% n = ^ ' 


(w-i)! 


(A-j)!^ 


vE<-W 


h=0 


\h 


Uh , 


' Wm - 


^(1) = 


(iV-j)! 


h=0 


Using them in 


«<,•>(/*) = 


u(h) 


i —1 


t=0 


Qi i) W-E(d B « ) (' i ) w<i ‘' ) ( ,i ) 


(ft = o, 1) 


gives equations (13.81) . (13.91) . where we denoted 

Qih := Re¬ 


using the above equations in (13.41) . we obtain the forms (j3.5j) and (13.6p for the coefficients 
Po, Pi, ■ ■ -,Pk- i and p m -i+ 1 , • • • ,p m -i,Pm, respectively. 

The remaining coefficients pi are to be determined so that 


m—l 


\Rr 


p m in = \\w -J2 p‘ b , 


12 
\L 2 


i=k 


W.= R n 


k—1 m \ 

E+ E U B r- 

J=0 j=m—l+ 1/ 


has the least value, where 


















WEIGHTED POLYNOMIAL APPROXIMATION OF RATIONAL BEZIER CURVES 


7 


Remembering that B™ and ]J ( . m " k ^ l> {k < i < m — l) are dual bases in the space II m’ l \ 
we obtain the formula 

Pi=(w,D\ mM } 

m—l / k—1 m \ 

j=k \j = 0 j = m —Z-f-l/ 


Observe that 


(R u , BJ') = Y,^r h ( B ',>. B T) 


h =0 
n 

= y^ j UJ h r h 

h =0 

n 

= y^ j UJ h r h 

h =o 


n \ / m\ [ m + n 


hj \j J \ h + j 


n \ (m \ [m + n 


-i 


<iAo JO 


V V i / V h + j 




Ih+j i 


where we use notation (13.101) . 

Using results of [20] and [12], we deduce that 


Kij :=(Bf, D[ m ’ k ’ l) } 

- f ( m V 1 (a + l + 1 )m-j ((3 + k + 1 )j 
\j )\i J ~ A:)!(m - l - i)\(a + l + 1 )i(/3 + k + l) k 

x di_ k (j — k; /3 + 2k, a + 21, m — k — l), 


where for h = 0,1,..., N we define 

(—l)' 1 TV! 


dh(x ; a, b, N ) = lim 


- V 

n^N n\(a + l)fc(6 + l)n-h “ 


i>=0 


{ — h)v { v + 1 — N) n _ v (b + l)n- v 
(n -\- o + b + 2) tv—d 
x Q n _ v (N - x; b, a + v + 1, TV - v — 1) 


with 


v 

Q P {t-,n,v,M) \= Y 

u =0 


(-p)u(p + p + p + l) u (-t) u 
u!(/i + 1) U (-M) U 


being Hahn polynomials (see, e.g., j9j §9.5]). It is easy to see that 


d h (x m , a, b, N ) 


(—x)jv + i(—iV - b) h 
(h - x)(a + l) h (6 + 1)tv 


(0 < h < N ). 
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Hence, some algebra gives 

T- = ( m \( m V 1 (~1 Y~ k (k - j)m-k-l +1 

VJ /I*/ {i - j){i - k)\{m - l - i)\ 

(o> + l + 1 ) m -j(f3 + k + l)j 
(a + l + 1 + k + l)j 

which is formula (13.111) . 

Now, (CUD readily follows. □ 


4. Implementation of the method 


4.1. Computing integrals (13.101) . Integrals (13.101) involving rational function cannot be 
evaluated exactly. However, we show that they can be computed numerically up to high 
precision using the method described in [ 8 j. 

Observe that formula (I3.10P can be written as 

(4.1) 4 = 2~ <t ~ n ( N \ J(a + N - l - h, 0 - k + h- V) (k < h < N - l; N := n + m), 


where 

(4.2) 


J(a, 6 ;i?) := / (1 — a:) a (l + x) b i!)(x) dx, 


'-i 


and 




(1 — x) l (l + x) k 


Notice that by assumption on the positivity of the weights Ui, the polynomial c j(t) has no 
roots in the interval [ 0 , 1 ], hence the function $(x) is analytic in a planar region containing 
the interval [—1, 1], This implies that the function $ can be well approximated by a sum 
of Chebyshev polynomials Tj(x) of the first kind: 


(4.3) 


Sm(x) = ItoJoOe) + 54 am 

3 = 1 


M 


(-1 < X < 1). 


See, e.g., [ 6 j. Clearly, 

J(a, b ; $) fa J(a, b ; Sm)- 

The coefficients 7 j are determined so that Sm interpolates i? at the abscissae ^ 
(0 < i < M ), hence 


cos(f7r/M) 


(4.4) 


7 3 '= 


2 — OjM 

M 


E"” 



cos 


Q7T 

M 


(0 <j<M). 


(The double prime on the sum means that the first and the last terms are to be halved.) 
The right-hand side of (14.41) is known to be efficiently computed by means of the FFT for 
real data (see [ 6 ], or [3], §5.1]; the authors recall that the FFT is not only fast, but also 
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resistant to roundoff errors). In practical implementation, the coefficients 7 j (0 < j < M ) 
are computed repeatedly for doubled values of M (M = 32, 64,...) until 

M 

( 4 - 5 ) M < 

i=M -3 

where £ is a prescribed tolerance. For better efficiency, in every step one should reuse all 
the previously computed values of the function d. 

Now, we have the following result. 


Algorithm 4.1 (Numerical computation of the integral J(a, 5; Sm ))• Given a, b > —1, let 
r \= b — a and s := a + b + 1. Define the sequence di (0 < i < M + 1) recusively by 


dM+i — dM 0 , 

(4-6) _ 2 r di + {i- s) d i+x - 2 7 * 

cti-i -:- 

1 + s 

Then we have 


(i = M,M- 1,...,1). 


(4.7) T(a, b ; Sm) — 2 s b B(a + 1 , b + 1 ) (70 — r do + s d x ). 

Proof. Function p{x) = (1 — x) a (l + x) b is a solution of the differential equation 

(1 — x 2 )p'(x) — [r — x(s — l)]p(x) = 0 , 
where r := b — a, s := a + b + 1. By p. Tlnn 2.1], we have 

/(l- x) a (l + x) b S M {x ) dx — (1 — x) a+1 (l + x) b+l T{x) + const (—1 < x < 1) 

for any continuous solution T of the differential equation 

(4.8) (1 — x 2 )T'{x) + [r — x(s + l)]T(x) = S M (x) (—1 < x < 1). 

Further, using the approach of m yields the recurrence relation 

(-i - s ) di_ 1 + 2 rdi + (i- s) d i+1 = 2 7* (i = 0,1,...) 

with d_i := d x and 7 j = 0 for j > M, satished by the Chebyshev coefficients di of the 
solution T of equation 04.81) . 

However, equation 04.81) has, in general, no continuous solution on [—1, 1]. On the other 
hand, by Pj. Tlnn 2.3] there always exists a polynomial solution of the modified equation 

(4.9) (1 - x 2 )T\x) + [r - x(s + l)}T(x) = S M (x) - \5 (-4 < x < 1), 

where 6 70 — r do + s d x . Obviously, 

J (a, b ; S M — = j (1 — x) a (l + x) 6 [ 5 'm(t) — |< 5 ] dx = 0, 

hence 


J (a, b■ S M ) = J (a, b ; |h) = \8 2 s B(a + 1,& + 1 ). 


□ 





10 STANISLAW LEWANOWICZ 1 ", PAWEL WOZNY, PAWEL KELLER 

4.2. Main algorithm. The presented method is summarized in the following algorithm. 


Algorithm 4.2 (Constrained polynomial approximation of the rational Bezier curve). 
Given the coefficients r 0 , ry,..., r n and u 0 , ay,..., uj n of the rational function 


R n (t ) 


Qn(t) 

u(t) 


Er=o Wi BN) 


(0 < t < 1), 


the coefficients po,Pi, ■ ■ ■ ,Pm of the polynomial 


P m (t) = Y J Pr B ?(t) (0<t<l) 

i=0 

minimising the error (I3.3[) with constraints (j3.4[) are computed in the following way. 


Step 1: Compute p 0 , pi,... ,pk- 1 by (13.51) . 

Step 2: Compute p m , p m - i, ■ ■ ■ ,Pm-i+ 1 by ([Ml). 

Step 3: Compute Cij(m, k, l, a, j3) for i,j — k,...,m — l by Algorithm ] 2. £1 
Step 4: Given £ > 0, compute the coefficients 7 j of the polynomial Sm (of. (14.31) ) by 
using FFT, with M determined so that (14.51) holds. 

Step 5: For h = k, k + 1,..., N — l, where N m + n, compute 

;= J(a + N — l - h, P — k + h\ S M ) 

by Algorithm ^.! . 

Step 6: For i — k, k + 1,..., m — l, compute p t by (13.7[) with I h replaced by 
/" := 2-°-» (Jfj J-f (k < h < N — l). 


5. Examples 


In this section, we present several examples of approximation of rational Bezier surfaces 
by Bezier curves with constraints, which we have described in Section [31 Computations 
were carried out on a computer with Intel Core i5 3.33GHz processor and 4GB of RAM. We 
used 18-digit arithmetic and set e := 10~ 16 in Step 4 of the Algorithm 14.21 to ensure that 
the integrals (I3.10P are computed within the accuracy close to the representation error. 

In the examples below, we use the notation 

E{t) := \\R n (t)-P m (t)\\ 

for the error function, and 

eoo := max E(t), e->(a, (3) ( [ (1 — t) a t^E 2 {t) dt 

* 6 [o, 1] ' Vi 0 

for the maximum and least-squares approximation error, respectively. 


Example 5.1 [Starling’s sketch] We consider the curve shown in Figure Q] (solid and 
red), obtained by joining two rational Bezier curves of degree eight each; we say that 
the composite curve has degree (8,8). The first curve is defined by the control points 
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(23,57), (-13,43), (29,58), (44,48), (30,42), (13,44), (-2,77), (42,83), (80,1), and the 
associated weights 1, 4, 3, 1, 5, 4, 7, 6, 1, while the second one by (80,1), (14,4), (3,54), 
(42,54), (51,42), (36,49), (66,12), (36,2), (47,3), and 1, 1, 4, 4,"2, 3, 3, 7, 8. The poly¬ 
nomial composite curve approximation of the above curve with a = j3 = |, constraints of 
C l - and C 2 -continuity, without and with subdivision, is shown in Figure [U □ 




Figure 1. Left: Bezier composite curve approximation (dashed and blue) of 
degree (13,8) to the rational composite Bezier curve (solid and red) of degree 
(8,8) with the end-point interpolation (k = l = 1). The errors are = 
{3.152,2.814} and e 2 (|, |) = {0.166,0.284}, respectively. Right: Bezier 
composite curve approximation (dashed and blue) of degree (12,11,7,6) to 
the same rational composite Bezier curve (solid and red) of degree (8,8), 
with one subdivision of each of its two parts, and the end-point derivative 
interpolation (k = / = 2). The errors are = {0.559, 0.811, 0.146, 0.231} 
and e 2 (- 7 , |) = {0.063, 0.104, 0.045, 0.081}, respectively. 


In Examples 15.21 and !5. 31 we compared our approach with the recently published methods 
of Huang et al. [7] and Lu HZ]. The idea of the first method is the following. Given the 
rational curve (13.11) . use degree elevation to obtain 

n-\-h 

J2^ h> rf h> Br h (t) 

= - (A =0,1,...), 

J2uf h> nr h (t) 

i =0 

and define the sequence of polynomial curves 

n-\-h 

tWt) := £ r< h > B ^ h (t) 
i =0 


{h = 0,1,...). 
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,p 

Then jU n+ h converges uniformly to as h —* oo, for any integer j > 0. The 

weakness of this approach is that the convergence may be rather slow. Also, notice the 
increasing degree of the approximating curves. 

In the iterative method of Lu HZ). the sequence of Bezier curves {V^} is constructed, 

n 

V n k (t)=J2 v < B - (*) ( ft = 0.1,...), 

i =0 

where 

v° = R n {ti) (i = 0,1,... ,n), 
with 0 = to < ti < ... < t n — 1, and 

Vi +l =Vi+ A (u? - V*(ti)) (i = 0,1,..., n; h = 0,1,...), 

A being a parameter. It is shown that 

lim V*(ti) = Rn(ti) {i = 0,1,..., n). 

h—too 

Also this process may be slow, even for carefully chosen factor A (cf. [16]). Another 
drawback of both methods is that only the simplest constraints (corresponding to k = l = 
1) are accepted. 

In the next two examples we let a = f3 = 0. 

Example 5.2 Our second test rational curve, which has the control points (14,1), (34, 25), 
(40,38), (—12,24), (5,21), (26,7), (18,41), (—13,34), (14,1) and the associated weights 
1, 3, 3, 4, 1, 7, 5, 3, 1, is shown in left part of Figure [2] (solid and red). We compared our 
approach with the methods of Huang et al. and Lu. We produced polynomial approxi¬ 
mation of degree m = 10 with end-point interpolation constraints (in our method we put 
k — l — 1). The results are shown in the left part of Figure [2] Table [2] below lists the 
errors of approximation in each case (with the number iter of iterations performed in the 
method of P3P- We see that the method of Huang et al, which is very simple and fast, 
gives much worse results than two other methods. The convergence of the method of Lu is 
slow: the result obtained after 100 iterations is 4 times less adequate than ours; also, the 
comparison of the execution times shows the advantage of our algorithm. □ 


Example 5.3 (17,12), (32,34), (-23,24), (33,62), (-23,15), (25,3), (30,-2), (-5,-8), 
(—5,15), (11, 8), and the associated weights 1, 2, 3, 6, 4, 5, 3, 4, 2, 1, is shown in the right 
part of Figure [2] (solid and red). We produced polynomial approximation of degree m = 10 
with end-point interpolation constraints. The results are shown in the right part of Fig¬ 
ure [2] and in Table [2] Again, we see that the methods of Huang et al. and Lu give much 
less adequate results than our method. □ 
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FIGURE 2. Polynomial approximation of degree 10 to the rational Bezier 
curve (solid and red) of degree 8 (left) and 9 (right), with end-point inter¬ 
polation (k = l = 1). Comparison of the present method (dashed and blue), 
Huang et al. method (dotted and brown) and Lu method (100th iterate; 
dotted-dashed and khaki). 

TABLE 2. Errors in the polynomial approximations to the rational Bezier 
curve, with end-point interpolation 



iter 

Example 15.21 

eoo e 2 (0,0) 

time[s] 

iter 

Example 15.31 

eoo e 2 (0,0) 

time[s\ 

Present method 

— 

0.664 

0.167 

0.15 

- 

0.398 

0.106 

0.15 

Huang et al. [7] 

— 

9.41 

3.98 

< 0.01 

- 

5.78 

3.03 

< 0.01 


25 

3.43 

0.850 

0.18 

25 

2.82 

0.653 

0.20 

Lu m 

50 

2.99 

0.729 

0.35 

50 

2.55 

0.584 

0.37 


75 

2.79 

0.671 

0.50 

75 

2.42 

0.550 

0.54 


100 

2.67 

0.636 

0.67 

100 

2.34 

0.529 

0.73 


Conclusions 

We present a method to solve the constrained least squares approximation of the rational 
Bezier curve by the Bezier curve. Important tools used are efficient evaluation of the Bezier 
coefficients of the constrained dual Bernstein basis polynomials associated with the Jacobi 
scalar product and numerical computation of some integrals involving rational functions. 
The new algorithm is particularly attractive when it is combined with the subdivision 
process. 
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